Astronomy & Astrophysics manuscript no. ccfbis 


©ESO 2011 


October 6, 2011 





Bisectors of the HARPS Cross-Correlation-Function* 
The dependence on stellar atmospheric parameters 

6. Bagturk 1 ' 2 , T. H. Dall 1 , R. Collet 3 , G. Lo Curto 1 , and S. O. Selam 4 



o 

(N 
-i— > 
O 

O 

in 



6 



> 
in 

O 



1 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Miinchen, Germany 

2 Ankara University, Astronomy & Space Sciences Research and Application Center, TR-06837 Ahlatlibel, Ankara, Turkey 

3 Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching bei Miinchen, Germany 

4 Ankara University, Faculty of Science, Department of Astronomy and Space Sciences, TR-06100 Tandogan-Ankara, Turkey 



Received date/ Accepted date 



ABSTRACT 



Context. Bisectors of the HARPS cross-correlation function (CCF) can discern between planetary radial-velocity (RV) signals and 
spurious RV signals from stellar magnetic activity variations. However, little is known about the effects of the stellar atmosphere on 
CCF bisectors or how these effects vary with spectral type and luminosity class. 

Aims. Here we investigate the variations in the shapes of HARPS CCF bisectors across the HR diagram in order to relate these to the 
basic stellar parameters, surface gravity and temperature. 

Methods. We use archive spectra of 67 well studied stars observed with HARPS and extract mean CCF bisectors. We derive previously 
defined bisector measures (BIS, v bot , c b ) and we define and derive a new measure called the CCF Bisector Span (CBS) from the 
minimum radius of curvature on direct fits to the CCF bisector. 

Results. We show that the bisector measures correlate differently, and non-linearly with logg and r cff . The resulting correlations 
allow for the estimation of log g and T e g from the bisector measures. We compare our results with 3D stellar atmosphere models and 
show that we can reproduce the shape of the CCF bisector for the Sun. 
Conclusions. 

Key words. Instrumentation: spectrographs - Techniques: radial velocities - Line: profiles - Stars: atmospheres - Stars: activity 
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1. Introduction 

Fundamental stellar parameters are best determined by study- 
ing high resolution, high signal-to-noise spectra. High resolu- 
tion spectroscopy is also essential for studying velocity fields in 
stellar atmospheres. There are a number of astrophysical phe- 
nomena leading to the emergence of macroscopic velocity fields 
in stellar atmospheres, which can be observed as spectral line 
asymmetries and radial velocity variations. Line bisectors, de- 
fined as the the loci of the midpoints on the horizontal lines ex- 
tending from one side to the other of spectral line profiles, are 
often employed to study the mechanisms that cause such asym- 
metries and variations. iGravl d2005l) related line asymmetries to 
a fundamental stellar parameter by showing that the blue-most 
point of the Fel /16252 line bisector is a very good indicator of lu- 
minosity class for the stars from late G to early K spectral types. 

Although it is possible to obtain high-resolution spectra with 
the use of state of the art telescope-spectrograph configurations, 
in order to achieve sufficient signal-to-noise ratio (SNR here- 
after), one typically has to either use long exposure times and/or 
combine many individual exposures. With this approach, how- 
ever, short-term line profile variations cannot be measured. To 
overcome this limitation, one can either average bisectors of in- 
dividual lines or employ cross-correlation functions (CCFs) in- 
stead. These functions are computed by cross-correlating ob- 
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servational stellar spectra with synthetic line masks with the 
objective of preserving line profile information; they can be 
thought of as near-optimal "average" spectral lines because the 
cross-correlation is based on lines of many different elements 
dDallet al.ll2006h . Because of this, CCFs have high SNR and 
represent one single instance of time. In addition, their compu- 
tation is faster compared to averaging bisectors of many individ- 
ual lines. On the other hand, averaging bisectors make it easy 
to select the individual lines according to line strength, excita- 
tion potential etc. However, one can also design cross-correlation 
masks by selecting lines with similar strength and/or excita- 
tion potential, or with a limited number of elements or ions. 
Designing masks other than the standard masks in the HARPS 
Data Reduction Software (DRS hereafter) is an ongoing study 
and will be the subject of another paper. 

Line bisectors delineate line asymmetries. In order to quan- 
tify the shape of bisectors, measures are defined. Bisector mea- 
sures frequently used in the literature are main ly based on the 
velocity span defined by Toner & Gravl (1 1988b as the velocity 
difference between two bisector points at certain flux levels of 
the cont inuum on a single lin e bisector, and the curvature de- 
fined bv lGrav & Hatzesl (1 19971) as the difference o f two velocity 
spans . Most of the extrasolar planet surv ey studies (lOueloz et al.l 
1 20011: iMartfnez Fiorenzano et al.l 120051) take the definition of 
iToner & Gravl (1988) as reference for the velocity span and 
adopt it to CCF bisectors using average velocities in two dif- 
ferent regions on the bisector, one at the top an d one close to 
the bo ttom instead of two real bisector points as IToner & Gravl 
(Il988l) did. They investigate if the variation of the bisector mea- 
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sure correlates well with that of the radial velocity. If such a 
correlation exists, then the radial velocity variation is attributed 
to stellar mechanisms instead of a planetary companion which 
would not cause variations in line asymmetry. 

It is lDall etalJ(l2006h that first gave a quantitative relation- 
ship between the stellar parameters (surface gravity, and abso- 
lute magnitude) and a combination of two bisector measures, 
namely the bisector inverse slope (BIS hereafter) and the curva- 
ture. Although their relation is based on a few data points, they 
stress that it can be possible to obtain the luminosity of a star 
from the shape of its average CCF bisector. Hence, the CCF bi- 
sector measures can be used as indicators of luminosity si milar 
to the blue-most point on classical single line bisectors as Grav 
d2005l) has shown. 

In this paper, we investigate the relationship between stellar 
parameters (specifically surface gravity and effective tempera- 
ture), and mean CCF bisector shapes of a number of late type 
stars observed with HARPS. In order to do this, we make use of 
the bisector measures frequently used in the literature as well as 
a new measure introduced here to represent the relative position 
of the upper segments of CCF bisectors. 

2. Observations and data analysis 

The da ta presented here a re obtained with the HARPS spectro- 
graph (M ayor et al.ll2003h which is installed on the 3.6 m tele- 
scope at the La Silla Observatory. HARPS provides high spectral 
resolution (1 15000 on average) and it is very efficient in provid- 
ing high radial velocity precision by keeping the instrumental 
profile (IP hereafter) stable thanks to its ultra-high internal sta- 
bility. 

We selected 67 late-type stars observed with HARPS (see 
Table [TJ, for which stellar parameters have been published pre- 
viously. In order to ensure a consistent set of stellar parameters 
across our data set, we preferably use references that determined 
these based on high-resolution spectra, and preferably obtained 
with HARPS . We cover the entire region of the HR diagram from 
late F to late K. Unfortunately, there are not many M stars in 
the HARPS archive with good SNR, and which have previously 
published trustworthy stellar parameters obtained by high reso- 
lution spectroscopy. Therefore, our study is limited to the spec- 
tral range from late F to late K stars. Our data set covers a wide 

range of surface gravities. 

Activity indices (Mount W ilson S-index (Vau ghan et alJ 

119781) and \ogR' HK dNoveslll984l) ) (see Table [Qi of all our pro- 
gram stars were determined with the DRS. 

As a rule of thumb, we eliminated spectra with SNR lower 
than 150 in the 60" 1 echelle order (~ 6025-6085 A) as calculated 
by the DRS. We inspected CCFs visually to guard against any 
obvious problems. HARPS CCFs are calculated for each of the 
72 echelle orders, which are then averaged and included by the 
DRS as a 73 rd order in the data products. We then extracted the 
bisectors of these averaged CCFs by taking the midpoints of hor- 
izontal lines connecting observation points on one wing to points 
on the other found by a cubic spline interpolation for each of the 
wings. We performed the procedure by using a new python code, 
developed by us. We obtained the radial velocities from the fits 
headers, that have been determined by the HARPS DRS from 
the centers of the Gaussian fits to the CCF profiles. 

In order to maximize the SNR, the mean bisector for each 
star has been computed, an example of which is given in Fig. Q] 
All the measurements presented in this paper have been per- 
formed on the mean bisector for each star. However, in order to 
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Fig. 1. Left: The CCF bisectors of HR 8658 (G0.5 V) computed 
using a G2 mask. Right: The mean CCF bisector. 
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Fig. 2. Mean bisectors of a Hor (K2 III) obtained with G2 and 
K5 masks. 



show the potential of our method in determining stellar parame- 
ters from a single high-resolution spectrum, we point out that the 
minimum SNR in the observations of HR 8658 (V - 6'" .64) is 
160, which has been reached in 230 seconds of integration time. 
For a 10 m .5 star, it is possible to reach an SNR of 150 in a one 
hour exposure. 

We used the CCFs computed by the HARPS pipeline (DRS). 
These functions are obtained by cross-correlating stellar spectra 
with a standard stellar mask matching as closely as possible the 
spectral type of the star. In Fig. [2] mean CCF bisectors of a Hor 
obtained by using the standard G2 and K5 masks of the HARPS 
DRS are given as an example. 

iDall et all d2006l) demonstrated that the shape of the CCF de- 
pends strongly on the mask used in the cross-correlation proce- 
dure (cf. their Fig. 3). In effect, the resultant CCFs obtained with 
different masks will be "averages" of different spectral lines. 
Thus, in order to take into account the effect of the mask when 
comparing CCF bisectors of different stars, we obtained mean 
CCF bisectors of all our program stars using both of the masks 
available in the HARPS DRS. This was accomplished during a 
data reduction mission to ESO Headquarters. 

2.1. Error Analysis 

The error on each bisect or point has been computed by following 
the procedure given bv lGravl(ll983h . The error in radial veloci- 
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Fig. 3. Left: The CCF bisectors of a Cen A (G2 V) computed 
using a G2 mask. Right: The mean CCF bisector. Note that the 
small scale structure of the bisector is constant and thus likely 
real. 

ties of our bisector points coming from photometric error can be 
expressed as, 



' V^SNRVT^CdF/dV);- 

where 1 / (SNR VI - term represents the photometric er- 
ror on flux for each bisector point i, where the line depth is indi- 
cated by ttj, and dF/dV is the slope of the profile at each point 
from which a bisector point is computed. The slope is obtained 
for each point from a linear fit to an array of points on the CCF 
profile, including the point itself and two adjacent points to it 
from both sides. 

The error on each bisector point is propagated on the bisector 
measures in a standard fashion by adding quadratically the errors 
of the bisector points used in the computation of the measure. 

2.2. Consistency of Bisectors 

Here we will discuss spectrum-to-spectrum variations of CCF 
bisectors, and their effects on the computation of mean bisectors 
and bisector measures. The RV for each spectrum has been cal- 
culated by the HARPS DRS by a Gaussian fit to the CCF profile. 
This RV has been subtracted from the individual CCF bisector. 
As a result, spectrum-to-spectrum variations in the CCF bisec- 
tor are mainly caused by (1) SNR differences due to differences 
in exposure time and variations in the conditions of the Earth's 
atmosphere, (2) instrumental effects, e.g. temperature drift, (3) 
errors coming from the reduction procedure, (4) variations in the 
stellar atmosphere due to magnetic activity, stellar oscillations, 
granulation, and any other variable velocity fields. Note that or- 
biting planets do not leave a signature on the CCF bisectors. 

In order to illustrate to what degree the bisectors are stable, 
we show the individual bisectors of a Cen A obtained during 
more than 10 hours of asteroseismological observations in the 
night of 21 April 2005 as well as their mean (Fig. [3). The SNR 
of the mean bisector is moderate (~875) compared to that for 
other stars in Table Q] all of which were observed for durations 
less than the total obervation time allocated for a Cen A in one 
night. 

To illustrate the long term consistency, 21 CCF bisectors 
of HD 21693 obtained in 21 different nights in the period 26 
December 2003 — 1 September 2009 is shown in Fig. |4] Fig. [3] 
and |4] show that the IP is very stable, both short term and long 
term. 
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Fig. 5. Left: Individual CCF bisectors of e Eri (K2 V) computed 
using a K5 mask. Right: The mean CCF bisectors for each night. 



An interesting spectrum-to-spectrum varia tion was observed 
for e Eri, which is a well-known active star (ICrolletalJ 120061 
Gray & Baliunas 1995). The star was observed on two nights, 
on 25 January 2007 and 2 February 2007. Five spectra were 
obtained on each of the nights. The change in individual CCF 
bisectors from one night to another is remarkable (Fig. |5), and 
illustrates both intra-night short term variation and longer term 
variations due to magnetic activity. In such cases, we used a 
mean CCF bisector obtained as an average of the bisectors ob- 
served in one night during which the shape is the most stable. 

2.3. Quantitative Bisector Measures 

A systematic change of CCF bisector shape with surface grav- 
ity and wi th effective te mperature is expected. This was demon- 
strated bv lGravl(l2005l) for classical single line bisector shapes 
over the HR diagram. In Sec. [3] we present the equivalent for 
the HARPS CCF bisectors. However, in order to quantify the 
variations, a measure that can numerically represent the bisector 
shap e is require d. 

lGravl (f2005). made use of the height of the blue-most point 
of a single line bisector as an empirical indicator of the luminos- 
ity of a star. However, CCF bisectors have shapes considerably 
different than those of single lines. The famous "C" shape of the 
bisector of the neutral iron line at /16253 in cool stars's spec- 
tra is physically attributed to granulation. Because granulation is 
depth dependent, the shape of the "C", and hence its blue-most 
point, changes from one star to another depending on how trans- 
parent the atmosphere is to granulation. Therefore, such a mea- 
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sure referring to the height of that point would be related to the 
surface gravity of the stars in our sample. However, as it is clear 
from the bisectors shapes given in both HR diagrams, there is not 
always a blue-most point on CCF bisectors, due to their nature as 
"average" spectral lines. In addition, if the HARPS IP is not sym- 
metric, this may affect the shap e of the bisector, and co uld cause 
a blue-most point to disappear. Living ston et al.l (1 19991) defined 
bisector amplitude as a proxy of the bisector shape to study line 
asymmetries observed in solar spectra. This measure also refers 
to an extreme blue point on a bisector. Furthermore, we will not 
be able to use the term flux deficit, which was defined and pro- 
po sed to be incr easing with effective temperature and luminosity 
bv iGravl J2010) because the physical context on which the def- 
inition is based, depends very strongly on absolute positions of 
sin gle spectra l lines. 

iDall et aD d2006l) made use of four bisector measures in their 
study on CCF bisectors. They adopted the definiton of veloc- 
ity span as the difference of average velocities between 10-40% 
and 55-90% of th e line depth from th e top, which is the bisec- 
tor inverse slope (lOueloz et al .11200 II) . leaving the term "veloc- 
ity span" for classical single line bisectors. They also adopted a 
curvature term, cj, which they define in segments, 20-30% (vi), 
40-55% (v 2 ), and 75-100% (v 3 ) of the lin e depth, as c b =((V3 - 
V2MV2-V1)), based on the definition of iPovich et a l. (2001). 
Additionally, they define the bisector bottom as the average ve- 
locity of the bottom four points on the bisector, minus the RV 
of the star. In order to represent the straight segment on the cen- 
tral part of CCF bisectors, they define another bisector measure 
which they call bisector slope. Because this measure is a variant 
of the BIS without any improvement, we will not use it. 

3. Results and Discussion 

In Fig. |6l mean bisectors of our program stars obtained by us- 
ing a G2 mask are illustrated on an HR diagram, based on their 
stellar parameters (Table [TJ depicted by asterisks on the half of 
the line depth. Mean bisectors on the HR diagram obtained by 
using a K5 mask are shown in Fig.|7]on the HR diagram with the 
same conventions as in Fig. [6] All measurements were obtained 
by using the standard HARPS G2 and K5 masks, respectively 

3.1. Correlations Between Mean Bisector Measures and 
Stellar Parameters 

A correlation between the bisector measures defined in Sec. I2.3I 
and surface gravity is obvious in Figs. [8] and [9] While BIS and 
Vbot both correlate well with surface gravity, cj is less convinc- 
ing. 

IDall et all d2006h gave a preliminary relationship between 
BIS+Cft and the surface gravity. Their relation is however based 
on only 9 solar type stars and biased by the effect of using both 
the G2 and the K5 masks. 

The bisector measures do not correlate with effective tem- 
perature as linearly as they do with surface gravity (Fig. [TOlfor 
bisectors obtained by using a G2 mask, and Fig. [TTJfor bisec- 
tors obtained by using a K5 mask), especially for the dwarf stars 
in the sample. The star, for which the bisector measures devi- 
ate most in the correlations, is Procyon. It is by far the hottest 
star in our sample. Unfortunately, we do not have intermediate 
temperature stars between Procyon and the next hottest star in 
the sample to trace the "non-linear" correlation hinted at in both 
Fig. [10] and inFig.QI] 

It is important to notice that the shape of the bisector also de- 
pends on the IP. Therefore, the correlations given here are valid 
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Fig. 6. Mean bisectors of program stars, computed by using a G2 
mask, on an HR diagram. Asterisks on the half line depth show 
the position of the star on the HR diagram based on its parame- 
ters obtained from the literatue and given in Table [TJ The scale 
for radial velocities spanned by mean CCF bisectors is given on 
the upper left corner. 



only for the HARPS instrument and cannot be used to estimate 
stellar parameters for other instruments. 
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Fig. 7. Mean bisectors of program stars, computed by using a K5 
mask, on an HR diagram. Same conventions apply as in Fig. [7] 
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Fig. 8. Surface gravity vs. bisector measures based on mean bi- Fig. 9. Surface gravity vs. bisector measures based on bisectors 
sectors obtained by using a G2 mask. Larger circle symbols in- obtained by using a K5 mask. Larger circle symbols indicate 
dicate larger effective temperatures. larger effective temperatures. 
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Fig. 10. Effective temperature vs. bisector measures obtained by Fig. 11. Effective temperature vs. bisector measures obtained by 
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Fig. 12. Mean CCF bisectors (plus signs) and appropriate NIST- 
Hahn fit functions (continuous curves) for the Sun, 27 Mon, HR 
2906, and a Mon from topleft in counter-clockwise direction. 



3.2. Direct fitting of the bisector and new quantitative 
measures 

In order to derive a bisector measure that is as objective as pos- 
sible, we have attempted to fit the bisector with a general ana- 
lytical function. We ha ve found tha t it is possible to make use of 
NIST-Hahn functions (lHahnl l l979) for fitting the entire bisector 
for a wide variety of shapes (Fig.[T2l. NIST-Hahn functions have 
the form 



f(x) = 



a + bx + cx 2 + rfx 3 
1 + fx + gx 2 + hx 3 



(2) 



where x is the flux, and f(x) is the velocity on the fit for a 
CCF bisector. Coefficients for the best fit for each mean CCF bi- 
sector has been obtained by Levenberg-Marquardt le ast-squares 
minim ization by using the python module MPFIT (Mar kwardl 

From the bisector shapes on the HR diagrams given in Figs. 
[6]and[7] it is expected that the relative position of the upper seg- 
ments of CCF mean bisectors changes systematically over the 
HR diagram. However, it is not possible to represent upper seg- 
ments of CCFs with a single velocity value as the bottom ve- 
locity represents the lower segment. Because the entire bisector 
can be fit by a NIST-Hahn function, the radial velocity value at 
a point on that function where the radius of curvature has a min- 
imum can be used as a measure of the relative position of the 
top of the bisectors. NIST-Hahn functions fit the overall bisec- 
tor shape but not irregularities on bisectors. This is advantageous 
especially for "wiggly" bisectors of F-type stars for which there 
is usually more than one point where the radius of curvature is 
low. Using NIST-Hahn functions; therefore, makes it possible to 
find a single point on the top of the bisectors where the radial 
velocity has a marginal value representing the top segment. 



Radius of curvature on a curve y, with one parameter x is 
given as 



R = 



>dx 2 > 



(3) 



We used Eq.[3]to analytically derive the radius of curvature at 
each point of the fit function (Eq. [2j) and find the velocity where 
the radius of curvature is minimum. We used this velocity value 
as a bisector measure and invesigated if this measure is corre- 
lated with stellar parameters. Fig. [13] and Fig. [14] show that this 
new bisector measure is a good indicator of the surface grav- 
ity. Similarly, the correlation between this measure and effective 
temperature is obvious in Fig. [15] and in Fig. [16] for both masks 
used in the cross-correlation procedure. Hence a positon of a star 
on the HR diagram is related to the radial velocity correspond- 
ing to the point where the radius of curvature of its bisector has 
a minimum. 

It is hinted in the bisector shapes in Fig. [6] and in Fig. [7] 
that, not only the positions of the upper and lower segments of 
mean CCF bisectors but also the radial velocity range spanned 
by them is related with the surface gravity and effective tempera- 
ture. Bisector In verse Slope (BIS) is a variant of the velocity span 
term defined by iToner & Gravl dl988h . Hence it is a measure of 
the velocity spanned by CCF bisectors in addition to its char- 
acteristic as a measure of the slope. However, BIS is based on 
the difference between average velocities between certain flux 
levels, which do not always represent the top and bottom seg- 
ments. Now that we have a good representation of the marginal 
velocities at the bottom with the bisector bottom and at the top 
with the RV at the minimum radius of curvature, we can define a 
good measure of velocity, that is spanned by mean CCF bisectors 
as the difference of these two measures. Since the ter m velocity 
span is reserved for classical line bisectors mostly (Pal l et alj 
2006), we name this measure as CCF Bisector Span (CBS). We 
do not make use of the radial velocities at the bottom segment 
of the NIST-Hahn function because the points at the bottom of 
bisectors are not well constrained due to low SNR and incompe- 
tency of the interpolation function. Instead, we use the average 
velocity of four of them (bisector bottom). Fig. [T3] [14] [15] and 
n~6l show that the CCF Bisector Span (CBS) is correlated very 
well with both surface gravity and effective temperature. 

3.3. Comparisons with Synthetic Spectra 

In order to better understand the properties and asymmetries of 
CCFs, we carried out a comparison of the observed CCF bisec- 
tors with theoretical ones which we extracted from synthetic 
spectra. More specifically, we computed synthetic solar spec- 
tra for the regions 5150 - 5200 A and 6215 - 6275 A using 
a realistic, three-dimensional (3D), time-dependent, radiative- 
hydrodynamic simulation of solar surface convection. The sim- 
ulation domain is a rectangular domain covering an area of 
about 6x6 Mm horizontally and extending for about 7 and 
6 pressure scale heights above and below the Sun's optical 
surface, respe ctively. The simulation was generated with the 
Stagger-Code (Nordlund & Galsgaard 1995), at a numerical res- 
olution of 240 3 . For the simulation, we employed a realis- 
tic equati on-of-state dMihalas et al.l Il988h . continuous opacity 
data from Gustafss on et al.l (1 19751) and Trampedach (2011, in 
prep.), and line opacity data from B. Plez (priv. comm.; see also 
iGustafsson et alj|2008l) . For the chem ical composition, we as- 
sumed the solar element mixture fromAspl und et aD {2009). 
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Fig. 13. Top: Surface gravity vs. radial velocity at the minimum 
radius of curvature. Bottom: Surface gravity vs. velocity span. 
Both measures are obtained on mean CCF bisectors which are 
computed by using a G2 mask. Larger circle symbols indicate 
larger effective temperatures. 
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Fig. 14. Top: Surface gravity vs. radial velocity at the minimum 
radius of curvature. Bottom: Surface gravity vs. velocity span. 
Both measures are obtained on mean CCF bisectors which are 
computed by using a K5 mask. Larger circle symbols indicate 
larger effective temperatures. 



We used the code SCATE dHaveket all 1201 lb to carry out 
a full spectral synthesis of the two spectral regions indicated 
above. For the calculations, we adopted a list of spectral lines 
with parameters extracted from VALD ( Piskunov et al. 1995; 



Rvabchik ovaetalJl997t iKupkaetal J 199^ 7 We selected a sam- 
ple of 20 solar simulation snapshots taken at regular intervals 
and covering in total one hour of simulated solar time. For each 
snapshot, we solved the monochromatic radiative transfer equa- 
tion under the approximation of local thermodynamic equilib- 
rium (LTE) for ^5 000 wavelength points in each spectral re- 
gion, for all grid points at the surface of the simulation, and 
along 16 inclined directions. Doppler shifts and broadening due 
to the presence of atmospheric velocity fields in the simulation 
are properly accounted for by the code, without the need to in- 



troduce additional micro- or macro-turbulence parameters. To 
model rotational broadening, however, we assumed a value of 
v sin i = 2 km/s for the Sun's rotational velocity. The final flux 
profiles were computed by averaging the spectra over time and 
over all grid-points at the surface of the simulation and by per- 
forming a disc-integration over the solid angle. Prior to the anal- 
ysis, the spectra were also rebinned and brought to the same 
wavelength resolution as the HARPS data. 

Our observed spectra of the Sun have been obtained from 
the daytime sky. We extracted the above mentioned wavelength 
regions from one of our high SNR spectra and normalized the 
fluxes in these wavelength intervals to the continuum level by 
using IRAF's standard packages. We then cross-correlated both 
the synthetic and observed solar spectra with two different masks 
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Fig. 15. Top: Effective temperature vs. radial velocity at the min- 
imum radius of curvature. Bottom: Effective temperature vs. ve- 
locity span. Both measures are obtained on mean CCF bisectors 
which are computed by using a G2 mask. Larger circle symbols 
indicate larger surface gravities. 
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Fig. 16. Top: Effective temperature vs. radial velocity at the min- 
imum radius of curvature. Bottom: Effective temperature vs. ve- 
locity span. Both measures are obtained on mean CCF bisectors 
which are computed by using a K5 mask. Larger circle symbols 
indicate larger surface gravities. 



("G" and "K"), specially designed for the same wavelength in- 
tervals. The masks available in the HARPS DRS pipeline are 
proprietary and not available outside the DRS, so they can only 
be used for the reduction of HARPS spectra. Hence, for these 
calculations, we decided to create new masks using VALD line 
lists for the spectral regions under study. For the sake of sim- 
plicity, we used the "extract stellar" function from VALD to re- 
trieve a list of lines with estimated central depths expressed as 
a fraction of continuum flux for two generic solar-metallicity 
stars belonging to the two spectral types. We constructed the 
masks by simply overlapping a series of impulse functions cen- 
tered at the wavelength of each spectral line and with amplitude 
proportional to the spectral line's depth. We were then able to 
cross-correlate both our synthetic and our observed spectra with 



the same masks. The cross-correlation was done using IDL's 
c .correlate function. Finally, we extracted the bisectors from 
the resultant CCFs of synthetic and observed spectra and com- 
pared them with one another. 

Fig. [T7] shows the CCF bisectors extracted from the syn- 
thetic (left panel) and observed (right panel) spectra in the 62 1 5- 
6275 A range. The CCFs were computed using the "G" mask, 
that is the more closely corresponding to the Sun's spectral type. 
The agreement between the CCF bisectors from synthetic and 
observed spectra is very good: the curvature and asymmetry of 
the observed bisector are well reproduced by the theoretical cal- 
culations. Fig. [18] shows the CCF bisectors from the synthetic 
and observed spectra for the same wavelength interval computed 
using the "K" mask instead. While the curvature of the CCF bi- 
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Fig. 17. Left: The bisector of the CCF obtained by the cross- 
correlation of the solar synthetic spectrum and "G" mask for the 
6215 - 6275 A region. Right: The bisector of the CCF obtained 
by the cross-correlation of a solar observed spectrum and "G" 
mask for the same region. Scale for x-axis is given on the figures 
and valid for both panels. 
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Fig. 19. Left: The bisector of the CCF obtained by the cross- 
correlation of the solar synthetic spectrum and "G" mask for 
5150 - 5200 A region. Right: The bisector of the CCF obtained 
by the cross-correlation of a solar observed spectrum and the 
same mask for the same region. The scale for x-axis is given on 
the figures and valid for both panels. 



1 .0 



0.8 



0.6 



0.4 







- 

50 m/s 





Radial Velocity (m/s) Radial Velocity (m/s) 

Fig. 18. Left: The bisector of the CCF obtained by the cross- 
correlation of the solar synthetic spectrum and "K" mask for the 
6215 - 6275 A region. Right: The bisector of the CCF obtained 
by the cross-correlation of a solar observed spectrum and "K" 
mask for the same region. Scale for x-axis is given on the figures 
and valid for both panels. 
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Fig. 20. Left: The bisector of the CCF obtained by the cross- 
correlation of the solar synthetic spectrum and "K" mask for the 
5150 - 5200 A region. Right: The bisector of the CCF obtained 
by the cross-correlation of a solar observed spectrum and the 
same mask for the same region. The scale for x-axis is given on 
the figures and valid for both panels. 



sector is more pronounced in the case of the synthetic spectrum, 
the overall agreement with the observed case is still good. Small 
differences of these kind are acceptable considering, e.g., that 
lines may be missing from the list used for the synthetic spec- 
trum calculations and that elemental abundances were not ad- 
justed to reproduce all line strengths. 

The corresponding CCF bisectors for the 5150 - 5200 A in- 
terval computed for the synthetic and observed spectra with the 
"G" and "K" masks are shown in Fig. [19] and [20] respectively. 

This region is characterized by the presence of one moder- 
ately strong and two strong Mg i lines at 5183, 5172, and 5167 A 
which carry a significant weight in determining the overall shape 
of CCF bisector in this region. The agreement between CCF bi- 
sectors derived from synthetic and observed spectra is actually 
excellent for this region, implying that the modelling of these 
strong lines in our theoretical calculations is robust and satisfac- 
tory. 



4. Conclusions 

We have shown that the well known CCF bisector parameters 
BIS, Vbot, and c;, correlate well with \ogg and T^. We have con- 
structed a new CCF bisector measure, the CBS, based on the 
point of least radius of curvature on the bisector, or alternatively, 
the point where the rate of change of the velocity fields is largest. 
The CBS measure too correlates very well with logg and T^. 
Although these correlations are not linear across the full range, 
we have shown that the shape of the CCF bisector changes in an 
understandable and predictable way across the HR diagram. 

We have shown that current state-of-the-art 3D stellar atmo- 
sphere models are able to reproduce the behavior of the CCF 
quite well for the Sun, which is very encouraging. Further mod- 
eling will be needed in order to extend this understanding to 
other late type stars across the HR diagram. 

The measures presented here have obviously predictive 
power; based on the CCF bisector, the fundamental atmospheric 
parameters of the star can be estimated. This may be relevant 
for studies of fainter stars where the SNR of individual spectral 
lines is not high enough for analysis. Such cases are still quite 
rare, partly because of the bias towards the study of brighter stars 
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with HARPS. Nevertheless, the results of this study can be use- 
ful in follow-up observations of exoplanet candidates discovered 
in transit surveys, by determining the parameters of especially 
the faint targets in a single observation based on the shape of its 
CCF bisector. 

So far, most planet searches have been trying to smear 
or average out the ef fects of the ste llar atmosphere (e.g., 
iDumu saue et al.l201 lb| a|). Only recently dBoisse et al.l201 lb has 
attempts been made to try to understand fully the effects of 
the stellar atmosphere with the aim of disentangling the effects. 
This is of particular importance in the context of the coming 
planet finding instruments, where sub-m/s precision will call for 
ways to discern and ultimately disentangle the effects of the 
stellar atmosphere from the planetary signal (e.g., ESPRESSO 
for the VLT and later Codex for the E-ELT: JPepe et al.ll20Tol 
iD'Odorico & The CODEX/ESPRESSO Teamll2007l) . 
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Table 1. P rogram stars. N is th e n umber of spectra ente ring the mean. mlAllende Prieto et al.1 d2004h; (2)ISousa et all d200l); (3)lGonzalez et al.1 
d2010h: (4llBruritt et alj f201Ch: C5)ICenarro et alJd2007h: (6)lda"Silva et alJfeOOd) : (7) ISanchez-Blazauez et al.1 d2006h : (8) lLuck & Heited d2007n : 
(9) ISantosetal.ld2004» : qOl lMc William! d!990h : (I D lHekker & Melendezl d2007f> : (12) lAbia et al.ldl988l) . 



HD 


Name 


Spec. 
Type 


N 


SNR 


Teff 

[K] 


[cgs] 


[Fe/H] 
dex 


S -index 


dex 


Ref. 




Sun 


G2V 


86 


850 


5777 


4.44 


0.00 


0.218 


-4.783 


1 


1461 


HR 72 


GOV 


38 


1500 


5765 


4.38 


0.19 


0.159 


-5.015 


2 


4915 


HD 4915 


GO 


32 


1850 


5638 


4.52 


-0.21 


0.209 


-4.782 


2 


7449 


HD 7449 


F9.5V 


29 


1800 


6024 


4.51 


-0.11 


0.178 


-4.780 


2 


10700 


rCet 


G8.5V 


438 


1275 


5310 


4.44 


-0.52 


0172 


-4.937 


2 


16417 


HR 772 


Gl V 


22 


2500 


5483 


4.16 


-0.13 


0.149 


-5.056 


2 


20794 


82Eri 


G8 V 


49 


1000 


5401 


4.40 


-0.40 


0.166 


-4.977 


2 


20807 


HR 1010 


GOV 


28 


1900 


5866 


4.52 


-0.23 


0.166 


-4.879 


2 


21019 


HR 1024 


G2V 


11 


2100 


5468 


3.93 


-0.45 


0.150 


-5.066 


2 


21693 


HD 21693 


G9IV-V 


21 


1050 


5430 


4.37 


0.00 


0.168 


-5.010 


2 


22049 


e Eri 


K2V 


12 


1000 


5133 


4.53 


-0.11 


0.527 


-4.500 


2 


23249 


<5Eri 


Kl III-IV 


12 


1000 


5150 


3.89 


0.13 


0.155 


-5.095 


2 


26967 


a Hor 


K2III 


150 


300 


4670 


2.75 


0.02 


0.108 


-5.457 


2 


38858 


HR 2007 


G4V 


22 


1400 


5733 


4.51 


-0.22 


0.173 


-4.925 


2 


38973 


HD 38973 


GOV 


15 


2400 


6016 


4.42 


0.05 


0.159 


-4.942 


2 


40307 


HD 40307 


K2.5V 


15 


650 


4977 


4.47 


-0.31 


0.198 


-4.996 


2 


50806 


HR 2576 


G5V 


14 


1150 


5633 


4.11 


0.03 


0.148 


-5.096 


2 


54810 


20 Mon 


K0III 


20 


1200 


4697 


2.35 


-0.30 


0.146 


-5.228 


2 


59468 


HD 59468 


G6.5V 


11 


1350 


5618 


4.39 


0.03 


0.162 


-5.004 


2 


60532 


HR 2906 


F6IV-V 


20 


4000 


6121 


3.73 


-0.25 


0.160 


-4.885 


3 


61421 


Procyon 


F5 IV-V 


70 


1550 


6485 


4.01 


0.01 


-0.170 




4 


61772 


HR 2959 


K3 III 


20 


1750 


3996 


1.47 


-0.11 


0.256 




5 


61935 


a Mon 


G9III 


20 


1350 


4879 


3.0 


-0.01 


0.134 


-5.266 


6 


65695 


27 Mon 


K2III 


5 


725 


4468 


2.3 


-0.14 


0.157 




6 


68978 


HD68978 A 


G0.5V 


25 


2200 


5965 


4.48 


0.04 


0.180 


-4.875 


2 


69830 


HR 3259 


G8V 


31 


1050 


5402 


4.40 


-0.06 


0.166 


-5. 008 


2 


72673 


HR 3384 


G9V 


12 


1150 


5243 


4.46 


-0.41 


0.180 


-4.930 


2 


73524 


HR 3421 


GOV 


10 


1500 


6017 


4.43 


0.16 


0.154 


-4.982 


2 


73898 


<TPyx 


G4III 


40 


950 


5030 


2.80 


-0.43 


0.115 


-5.294 


7 


81169 


A Pyx 


G7 III 


5 


1200 


5119 


3.13 


-0.07 


0.142 


-5.148 


8 


81797 


a Hya 


K2 II-III 


5 


400 


4186 


1.7 


0.00 


0.212 




6 


84117 


HR 3862 


F8 V 


25 


1750 


6167 


4.35 


-0.03 


0.178 


-4.824 


9 


85444 


39 Hya 


G6-G8III 


5 


3000 


5000 


2.93 


-0.14 


0.245 


-4.853 


10 


85512 


HD 85512 


K6V 


6 


1450 


4715 


4.39 


-0.32 


0.399 


4.858 


2 


85859 


HR 3919 


K2III 


10 


575 


4340 


2.35 


-0.03 


0.150 




10 


90156 


HD 90156 


G5V 


4 


1250 


5599 


4.48 


-0.24 


0.170 


-4.952 


2 


96700 


HR 4328 


GOV 


8 


1650 


5845 


4.39 


-0.18 


0.164 


-4.939 


2 


98430 


<5Crt 


K0III 


15 


800 


4580 


2.35 


-0.43 


0.145 




11 


100407 


£Hya 


G7III 


30 


1500 


5045 


2.76 


0.21 


0.212 


-4.947 


4 


102365 


HR 4523 


G2V 


60 


1100 


5629 


4.44 


-0.29 


0.187 


-4.792 


2 


102438 


HR 4525 


G6V 


37 


1400 


5560 


4.41 


-0.29 


0.170 


-4.951 


2 


102870 


/3 Vir 


F8V 


28 


1850 


6050 


4.01 


0.12 


0.181 


-4.823 


4 


114613 


HR 4979 


G3V 


75 


850 


5729 


3.97 


0.19 


0.165 


-4.949 


2 


115617 


61 Vir 


G7V 


36 


900 


5558 


4.36 


-0.02 


0.167 


-4.984 


2 


115659 


y Hya 


G8III 


15 


950 


5110 


2.90 


0.03 


0.133 


-5.203 


11 


123123 


jt Hya 


Kl III-IV 


21 


725 


4670 


2.65 


-0.16 


0.141 


-5.355 


11 


128620 


a Cen A 


G2V 


982 


875 


5745 


4.28 


0.22 


0.160 


-4.948 


4 


128621 


aCenB 


K1IV 


6 


825 


5145 


4.45 


0.30 


0.232 


-4.820 


4 


130952 


11 Lib 


G9III 


30 


2000 


4820 


2.91 


-0.38 


0.148 


-5.180 


10 


134060 


HR 5632 


GOV 


44 


2050 


5966 


4.43 


0.14 


0.158 


-4.980 


2 


134606 


HD 134606 


G6IV 


8 


1500 


5633 


4.38 


0.27 


0.150 


-5.070 


2 


134352 


HR 5699 


G4V 


58 


2000 


5664 


4.39 


-0.34 


0.167 


-4.968 


2 


136713 


HD 136713 


K2V 


28 


750 


4994 


4.45 


0.07 


0.311 


-4.750 


2 


144585 


HR 5996 


G1.5V 


17 


1500 


5914 


4.35 


0.33 


0.148 


-5.042 


2 


147675 


y Aps 


G8III 


25 


1800 


5050 


3.50 


0.05 


0.248 


-4.834 


12 


154962 


HR 6372 


G6IV-V 


5 


2300 


5827 


4.17 


0.32 


0.139 


-5.158 


2 


160691 


fj. Ara 


G3 IV-V 


279 


400 


5780 


4.27 


0.30 


0.149 


-5.101 


2 


162396 


HR 6649 


F9V 


23 


2000 


6090 


4.27 


-0.35 


0.152 


-4.992 


2 
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Table 1. Pr ogram stars. N is th e n umber of spectra ente ring the mean. (l)| Allende Prieto et al. (2004); (2) Sousa et al. (2008); (3) Gonzale z et al.l 
( 120101); (4)lBruntt et al] f20Toh; f5)ICenarro et alj d2007l); (6)lda Silva et all J2006t) : (7) ISanchez-Blazauez et al.l d20o"6h ; ( 8) lLuck & HeiteJ d2007n : 
(9) ISantosetal.ld2004l) ; (lO) lMcWilliaml d 19901) ; (1 D lHekker & Melendezl ( 120071) ; (T2) lAbia et al.ldl988l) . 



HD 


Name 


Spec. 


N 


SNR 


T e ff 




[Fe/H] 


S -index 




Ref. 






Type 






[K] 


[cgs] 


dex 




dex 




i / ivyu 


rlK OW 


ro V 


1 1 
1 1 


jUUU 


AH/1 < 


A 1 A 


U.UO 


U. 14Z 


-J.U / o 


L 


176986 


HD 176986 


K2.5V 


28 


600 


5018 


4.45 


0.00 


0.264 




2 


179949 


HR 7291 


F8V 


8 


2400 


6287 


4.54 


0.21 


0.192 


-4.747 


2 


190248 


£Pav 


G8IV 


17 


1300 


5604 


4.26 


0.33 


0.152 


-5.070 


2 


207129 


HR 8323 


G2V 


42 


1050 


5937 


4.49 


0.00 


0.199 


-4.795 


2 


209100 


elnd 


K5V 


6 


500 


4754 


4.45 


-0.20 


0.419 


-4.767 


2 


210918 


HR 8477 


G2V 


12 


1600 


5755 


4.35 


-0.09 


0.160 


-4.985 


2 


215456 


HR 8658 


G0.5V 


29 


2000 


5789 


4.10 


-0.09 


0.146 


-5.068 


2 


221420 


HR 8935 


G2IV-V 


9 


2150 


5847 


4.03 


0.33 


0.135 


-5.188 


2 



